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^ l' We present low temperature tunneling density of states measurements of Pauli-limited Al films 

D ' in which the Zeeman and orbital contributions to the critical field are comparable. We show that 

, films in the thickness range of 6-7 nm exhibit a reentrant parallel critical field transition which is 

. associated with a high entropy superconducting phase, similar to the high entropy solid phase of 

' ''He responsible for the Pomeranchuk effect. This phase is characterized by an excess of states near 

I 1 1 the Fermi energy so long as the parallel critical field transition remains second order. Theoretical 

P3 , fits to the zero bias tunneling conductance are in good agreement with the data well below the 

O ' transition but theory deviates significantly near the transition. The discrepancy is a consequence of 

O ' the emergence of e-e interaction correlations as one enters the normal state. 

U ■ 

p PACS numbers: 74.25.Dw,74.78.Db,74.50.+r 
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I. INTRODUCTION 



^ ' Phase transitions are ubiquitous in physics and continue to be a large and important part of current condensed 
'"^ matter research. Perhaps the most common example of a classical phase transition is one in which a system undergoes a 
^ first-order transition from a disordered phase to a phase exhibiting long range order, such as a hquid to solid transition. 
O Because liquid has a higher entropy, latent heat must be removed from the system before the solid phase can form. 
I Though examples are exceedingly rare in nature, systems which undergo a transition from a disordered phase to an 

ordered phase of higher entropy exhibit a number of counterintuitive and fascinating properties such as reentrance and 
^-H , negative latent heats. The most famous example is the liquid to solid transition in ■^He. Subtle correlations between its 
^ ' spin and coordinate degrees of freedom produce an anomalously low entropy in liquid ■^Hc. Consequently, the entropy 
, of solid '^He is higher than that of the liquicii near the melting curve. As was first pointed out by Pomeranchuk;^ the 
' resulting negative latent heat can be used to produce a cooling effect if one converts liquid to solid by increasing the 
Qv^ \ pressure. This is commonly known as the Pomeranchuk effect and it was, in fact, the method used to cool liquid ■^He 
■ in the series of experiments that led to the discovery of superfluid '^Heii The only other quantum system known to 
' exhibit a high entropy ordered phase is Pauli-limited superconductivity."^! In analogy to "^He, the anomalous entropy 
is associated with the interplay between electronic spin and orbital degrees of freedom where the superconducting 
phase plays the role of the solid, the normal phase the liquid, and parallel magnetic field the role of pressure. In the 
present paper, we address the possibility of realizing a superconducting version of the Pomeranchuk effect. 

Magnetic fields generally have a detrimental effect on superconductors via two independent channels. The first 
is an orbital effect associated with the fact that cyclotron motion is incompatible with the formation of Cooper 
pairs and hence superconductivity; for the vast majority of superconducting systems the critical field transition is 
Q I completely dominated by the orbital response of the conduction electrons. If a magnetic field is applied in the plane 
O • of a superconducting film whose thickness is much less than the superconducting coherence length and whose electron 
diffusivity is low,'^ then the critical field transition is, in fact, mediated by the Pauli spin polarization of the electrons: 
the system will undergo a well defined first-order phase transition to the normal state when the electron Zeeman 
' splitting is of the order of the superconducting gapi^ When realized in low atomic mass superconductors such as 
M I Al and Be, this "spin-paramagnetic" (S-P) transition, as it is commonly known, exhibits a variety of novel phases 



along with some very unusual dynamics ^(iSiSiiS Here we present tunneling measurements of the density of states 
(DoS) in Al films which are about 2-3 times thicker than those typically used in S-P studiesiii In this marginally 
Pauli-limited regime, the critical field behavior of the films can become reentrantii^ The reentrance occurs because 
the superconducting state has higher entropy than the normal state. This excess entropy is expected to be reflected 
in an enhancement of the Fermi energy DoS above that of the normal state .^d^ds 

In a classic S-P geometry a planar magnetic field is applied to a film of thickness much less than the Pippard 
coherence length ^. The relative importance of the orbital response compared with the spin polarization can be 
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quantified by the dimensionless orbital depairing parameter— 

D{ed)^A 



6hii 



where D is the electron difFusivity, d is the film thickness, e is the electron charge, Aq the zero-temperature, zero-field 
gap energy and /i b the Bohr magneton. Fulde and coworkers^ showed that in the absence of spin-orbit coupling, the 
low temperature critical field behavior would remain first-order for c < 1.66, but interestingly, for (thicker) films with 
values in the range 1.66 < c < 2.86, the critical field behavior was predicted to be reentrant. Reentrance has indeed 
been observed in Ali4 and Ti]>ii^ films of comparable thickness to those used in this study. Tunneling DoS studies, 
however, have been quite limited and though considerable theoretical work has been done on the problemji2ii& to the 
best of our knowledge a quantitative comparison between theory and experiment was never made in the regime where 
the DoS enhancement is present. 



II. EXPERIMENTAL DETAILS 



A series of Al films with thicknesses ranging from 6 to 8 nm and corresponding 100 mK normal state sheet resistances 
ranging from i? = 10 to 20 were fabricated by e-beam deposition of 99.999% Al onto fire polished glass microscope 
slides held at 84 KiiS The depositions were made at a rate of ~ 0.2 nm/s in a typical vacuum P < 3 x 10^^ Torr. 
After deposition, the films were exposed to the atmosphere for 0.5-4 hours in order to allow a thin native oxide 
layer to form. Then a 14-nm thick Al counterelectrode was deposited onto the film with the oxide serving as the 
tunneling barrier. The counterelectrode had a parallel critical field of ~ 1.1 T due to its relatively large thickness, 
which is to be compared with H^w ~ 3 T for the films. The junction area was about 1 mmxl mm, while the junction 
resistance ranged from 10-20 kfi depending on exposure time and other factors. Only junctions with resistances much 
greater than that of the films were used. Measurements of resistance and tunneling were carried out on an Oxford 
dilution refrigerator using a standard ac four-probe technique. Magnetic fields of up to 9 T were applied using a 
superconducting solenoid. A mechanical rotator was employed to orient the sample in situ with a precision of ~ 0.1°. 



III. RESULTS AND DISCUSSION 



Shown in the upper panel of Fig. ^ is the temperature dependence of the parallel critical field of a 2.5 nm thick Al 
film having Tc — 2.7 K, Ao/e = 0.41 mV, ^ 15 nm, and sheet resistance R ^ 1 kfi. Using the bulk DoS for Al 
i^o = 2 X lO^'^ (eVcm'^) to estimate D in Eq. J^l, the orbital depairing parameter for this film is c = O.Q2^ The 
upper panel of Fig.^represents the classic S-P phase diagram in which a high temperature line of second order phase 
transitions terminates into a line of first order transitions at the tricritical point Ttri ^ 0.3Tcii^ The low temperature 
superconducting (S) and normal (N) phases are separated by a robust coexistence region in which the state of the 
system is solely determined by the system state prior to entering the region, i.e. the state memory region (SM)j2i In 
the lower panel of Fig. ^ we show the corresponding phase diagram for a 6 nm thick Al film having Tc = 2.1 K, 
Ao/e = 0.32 mV, f = 56 nm, and i? = 20 $7. Note that this thicker film has a substantially lower critical field as a 
consequence of orbital depairing, and that the structure of its phase diagram is dramatically altered. The critical field 
behavior is now reversibly reentrant by virtue of the local maximum in the parallel critical field near T/Tc = 0.45. 
Indeed, a naive estimate of the orbital depairing parameter c ~ 2.2 falls within the reentrant range. In the upper 
inset we show that the resistive parallel critical field at 400 mK is higher than it is at 60 mK. Reentrance is easily 
demonstrated by cutting across the phase diagram along the dashed line in Fig. ^ as is shown in the lower inset. 
The numbers in this curve represent the sequence of temperature scans at a constant field. The reentrant region was 
entered at 60 mK by lowering the field from a supercritical value of 3.5 T to the field at which the temperature scans 
were performed, 3.320 T. The low temperature hysteresis is a consequence of the fact that the transition becomes 
first order below ^ 0.4 K . 

Shown in Fig.|21are tunneling spectra at several values of parallel critical field. At the temperatures at which these 
data were taken the tunneling conductance is simply proportional to the electronic DoSi^Si The Zeeman splitting of 
the BCS coherence peaks can be seen in the 2.2 T curveiSS As the field is increased, the gap begins to collapse due 
to orbital pair breaking, and near the critical field transition the subgap Zeeman features reach the Fermi energy, 
thereby producing a zero bias peak in the DoS curve. This superconducting DoS enhancement is evident in the 3.02 
T curve of Fig. |21 where the DoS at V = is clearly higher than that of the 3.06 T normal state curve. 

In Fig. |3] we plot the zero-bias tunneling conductance as a function of parallel field at a variety of temperatures 
for a 7 nm Al film of resistance R = 15.8 f2. The tunneling conductance has been normalized by the normal 
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state value. The arrows indicate the onset of superconductivity. Because of the reentrance effect, the onset is not 
monotonic in temperature. Suzuki and coworker a^^i^^ reported a 1% DoS enhancement over only a very narrow 
range of temperatures 550 mK < T < 650 mK in Al films in the thickness range of 7-9 nm. In contrast to these 
earlier reports, the DoS enhancement peak appears below 500 mK and, in fact, grows in magnitude with decreasing 
temperature. This film did not exhibit a tricritical point, therefore we believe that the transition remained second 
order down to 60 mK. Theoretically one can compute such curves, assuming that the film's parameters are known, by 
numerically solving the self-consistent equations for the order parameter and the "molecular" magnetic field together 
with the Usadel equations for the semiclassical Green's functions; these equations are derived in Ref . Il3l and we refer 
the reader to this work for the details (see also Ref. for numerical calculations of the temperature dependence 
of the critical field). Here we remind that the input parameters needed to solve these mean- field equations at a 
given temperature T are: the orbital pair-breaking parameter c (cf. Eq. the spin-orbit scattering parameter 

b = h/3AoTso, where Tso is the spin-orbit scattering time; the Fermi-liquid parameter G'^ which, in the normal state, 
describes the renormalization of the spin susceptibility; and the gap Aq. In Fig.^jwe compare our data to theoretical 
curves calculated using parameters values c = 0.79, b — 0.052, G" = 0.155, and the measured gap Ao = 0.32 meV, 
for three different temperatures. We note that the value of b is in agreement with what is found in the literature,^ 
while G'^ is consistent with what was previously measured in experiments performed in the normal state on similar 
samplesi^ On the other hand, the value of c used to fit the data is significantly lower than that estimated by Eq. 
and ostensibly below the reentrant threshold discussed above. However, the range of c over which reentrance can be 
observed is affected by both spin-orbit scattering and the Fermi-liquid effect, as discussed in detail e.g. in Refs. .SlISI 
Furthermore, Eq. (QJ is valid only in the "local" limit £ d, where £ is the mean free path for impurity scattering. 
A more general relation for c is found by multiplying Eq. J^l by a function f{£/d), which describes the cross-over 
between the "local" and "non-local" {£/d ^ 1) electrodynamics in thin films, ^— and whose asymptotic behavior is 
f{x) ~ 3/4a; for a; ^ 1; from the measured conductance of our sample, the bulk DoS z/q, and the bulk Fermi energy 
Ep = 11.7eV, we estimate that f/d « 3.4 and hence c ~ 0.73, in good agreement with the numerical fitting; therefore 
the values of all the three fitting parameters c, b and G° are consistent with independent measurements or estimates. 

At fields well below the critical field there is excellent agreement between theory and experiment. As the field 
increases, however, the measured zero-bias tunneling DoS falls below the computed value by an amount we define as 
r] in Fig. 21 The origin of rj cannot be attributed to superconducting fluctuations near the critical field, as the Cooper 
pairing contributions to the DoS are shifted by the parallel magnetic field to energies of the order of the Zeeman 
energy;^ it can, however, be explained as follows: increasing the field depresses the order parameter, so that more 
electrons can contribute to interaction correlations which are not included in the mean- field analysis of Ref. ^3 these 
correlations cause a renormalization of the DoS close to the transition. In fact, it is well known that e-e interactions 
lead to a logarithmic zero bias correction to the 2D tunneling DoS, Siy oc In(TT), where r is the elastic scattering 
timem^ In the inset of Fig. 0] we plot ij as a function of InT. Note that above 300 mK, rj is well described by a In T 
temperature dependence. The low temperature saturation is partly a consequence of the ~ 10 probe voltage used 
to measure the tunnel conductanceiS^ The horizontal dashed line in the inset represents the magnitude of the ZBA 
as obtained by a direct measurement of the normal state DoS at 60mK. Note that the ZBA and ry are of the same 
magnitude and that a close inspection of the data in Fig. ^ reveals that the region over which there is disagreement 
between mean-field theory and experiment grows with increasing temperature, as would be expected. 

In Fig. [S] we plot the magnitude of the DoS peak as a function of temperature for a 7 nm and a 6 nm film. 
Interestingly, the temperature dependence of the peak is quite different for the two films. The critical field behavior 
of the 7 nm film did not show any signs of hysteresis down to 60 mK, indicating that the transition remained second 
order. In contrast, the 6 nm film exhibited a tricritical point at T/Tc ~ 0.2, as can be seen in Fig. ^ Note that 
the DoS peak in this film has a local maximum at the tricritical point, and that the DoS enhancement is rapidly 
extinguished as the temperature is lowered below Ttri- This may explain why earlier reports of the DoS enhancements 
by Suzuki and coworkers^ were limited to a relatively narrow range of temperatures, i.e. temperatures lying between 
a tricritical point and the maximum of the critical field curve. 

In summary, we show not only that marginally thick Al films exhibit reentrance, but also that the associated DoS 
enhancement can be observed down to mK temperatures so long as the critical field transition remains second order. 
Since the Pomeranchuk effect is based on the latent heat of a first order transition, it seems unlikely that one can 
produce cooling in Al films. Nevertheless, the fact that the DoS enhancement peak occurs at the Fermi energy affords 
one an unparalleled opportunity to study the interplay between Coulombic e-e interaction fluctuations and pairing 
fluctuations in the vicinity of the transition. Indeed, it is fortuitous that the suppression of the DoS associated with 
the former is comparable in magnitude to the enhancement of the latter. 
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FIG. 1: Upper panel: phase diagram of a 2.5 nm thick Al film (Tc = 2.7 K) in parallel field, where S is the superconducting 
phase, N the normal state, SM is the state memory coexistence region, and Tu-i is the tricritical point. Lower panel: phase 
diagram of a 6 nm thick Al film {Tc = 1.9 K). Note the local maximum in the second-order portion of the critical field 
curve. The horizontal dashed line represents a reentrant cut through the phase diagram. Upper inset: Two resistive critical 
field transitions in the 6 nm film showing a higher Hc\\ at higher temperature. Lower inset: Reentrant superconducting 
transition corresponding to the cut through the 6 nm phase diagram represented by the dashed line. The numbers represent 
the temperature scan sequence. 
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FIG. 2: Tunneling density of states spectra of a 7 nm Al film in several parallel magnetic fields. Note that the zero bias peak 
in the 3.02 T curve rises above the 3.06 T normal state curve. The modest suppression of the DoS at zero bias in the normal 
state curve is an e-e interaction effect and is commonly known as the zero bias anomaly. 
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FIG. 3: Normalized density of states at zero bias for a 7 nm thick film. The peaks correspond to excess states at the Fermi 
energy associated with the onset of superconductivity. Note that the onset critical fields, indicated by the arrows, are not 
monotonic in T. 
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FIG. 4: The dashed lines represent theoretical fits to the zero bias tunneling data at 60 mK, 300 mK, and 500 mK, top to 
bottom. The respective curves have been shifted for clarity. Only that portion of the data well below the critical field was used 
in the fits. The magnitude of the discrepancy between the theory and data near the transition is represented by the parameter 
rj. Inset: magnitude of the offset between theory and data in the normal state as a function of temperature for the 7 nm Al 
film. The solid line is provided as a guide to the eye. The horizontal dashed line represents the relative magnitude of the zero 
bias anomaly in the film. 
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FIG. 5: Magnitude of DoS enhancement peaks such as those in Fig. |21as a function of reduced temperature for a 6 nm and 
a 7 nm thick AI film. The 7 nm film did not exhibit a finite tricritical point. The 6 nm film critical field transition became 
first-order below T/T^ ~ 0.2. 



